Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 30 November 2009 (MN WI^ style file v2.2) 



A comprehensive population synthesis study of 
post-common envelope binaries 



o 
o 

> 

o 

o 

m 

^' 

o 



(N 
> 
(N 



m 
o 

o 



p. J. Davis,! U. Kolb,i B. Willems^ 

-"^ The Open University, Department of Physics and Astronomy, Walton Hall, Milton Keynes MK7 6AA 
'^Northwestern University, Department of Physics and Astronomy, 2131 Tech Drive, Evanston, IL 60208, USA 



30 November 2009 



ABSTRACT 

We apply population synthesis techniques to calculate the present day population of 
post-common envelope binaries (PCEBs) for a range of theoretical models describing 
the common envelope (CE) phase. Adopting the canonical energy budget approach we 
consider models where the ejection efficiency, acE is either a constant, or a function 
of the secondary mass. We obtain the envelope binding energy from detailed stellar 
models of the progenitor primary, with and without the thermal and ionization energy, 
but we also test a commonly used analytical scaling. We also employ the alternative 
angular momentum budget approach, known as the 7-algorithm. We find that a con- 
stant, global value of acE ^0.1 can adequately account for the observed population 
of PCEBs with late spectral-type secondaries. However, this prescription fails to re- 
produce IK Pegasi, which has a secondary with spectral type A8. We can account for 
IK Pegasi if we include thermal and ionization energy of the giant's envelope, or if 
we use the 7-algorithm. However, the 7-algorithm predicts local space densities that 
are 1 to 2 orders of magnitude greater than estimates from observations. In contrast, 
the canonical energy budget prescription with an initial mass ratio distribution that 
favours unequal initial mass ratios {n{qi) oc q^^'^^) gives a local space density which is 
in good agreement with observations, and best reproduces the observed distribution 
of PCEBs. Finally, all models fail to reproduce the sharp decline for orbital periods, 
-forb ^ 1 d in the orbital period distribution of observed PCEBs, even if we take into 
account selection effects against systems with long orbital periods and early spectral- 
type secondaries. 
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1 INTRODUCTION 



The common e nvelope (CE) phase was proposed by 
iPaczvnskil l| 19761 ) to explain the small orbital separation of 
compact binaries such as cataclysmic variables (CVs) and 
double white dwarf binaries, because a much larger orbital 
separation of between approximately 10 and 1000 R© was 
required in order to accommodate the giant progenitor of 
the white dwarf. (For reviews see Iben & Livio 1993, Taam 
& Sandquist 2000). 

In the pre-CE phase of evolution the initially more mas- 
sive stellar component (which we henceforth denote as the 
primary) evolves off the main sequence first. Depending on 
the orbital separation of the binary, the primary will fill its 
Roche lobe on either the giant or asymptotic giant branch 
and initiate mass transfer. If the giant primary possesses a 
deep convective envelope (i.e. the convective envelope has a 
mass of more than approximately 50 per cent of the giant's 
mass; Hjellming & Webbink 1987) the giant will expand in 



response to rapid mass loss. As a result, the giant's radius 
expands relative to its Roche lobe radius increasing the mass 
transfer rate. As a consequence of this run-away situation, 
mass transfer commences on a dynamical timescale. The 
companion main sequence star (henceforth the secondary) 
cannot incorporate this material into its structure quickly 
enough and therefore expands to fill its own Roche lobe. 
The envelope eventually engulfs both the core of the pri- 
mary and the main sequence secondary. 

The friction the stellar components experience due to 
their motion within the CE extracts energy and orbital an- 
gular momentum from the orbit causing the orbital separa- 
tion to decrease. If enough of the orbital energy is imparted 
on to the CE before the stellar components merge, then the 
CE may be ejected from the system leaving the core of the 
primary (now the white dwarf) and the secondary star at a 
greatly reduced separation. 

Despite extensive three dimensional hydrodynamical 
simulations (e.g. Sandquist, Taam & Burkert 2000), the 
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physics of the CE phase remains poorly understood, includ- 
ing the efficiency with which the CE is ejected from the 
system. Due to the difficulty in modelling the CE phase, bi- 
nary population synthesis calculations commonly resort to 
describing the CE phase in terms of a simple energy budget 
argument. A fraction acE of the orbital energy released as 
the binary tightens, Ai?orb, is available to unbind the giant's 
envelope from the core. Hence if the binding energy is -Bbind 
we have 



Sbin 



aCEAiSorb 



(1) 



(e.g. de Kool 1992; de Kool & Ritter 1993; Willems & Kolb 
2004). The efficiency qce is a free parameter with values 
< acE < 1, albeit values larger than unity are discussed, 
dep ending on what t e rms a re included in the quantity -Ebind- 
iNelemans et all (|2000l) used this approach to recon- 
struct the common envelope phase for the observed sample 
of double white dwarf binaries. They found that acE < 
for the first phase of mass transfer, which is unphysical. This 
led them to test an alternative description of the CE phase, 
based on the change in the orbital angular momentum of the 
binary, A J, during the C E phase. Their '7-algorithm' foll ows 
the prescription given bv IPaczvriski fc ZiolkowskH (|l967h . 



AJ = 7 



11" 

Mb 



(2) 



where Mcnv is the mass of the giant's envelope, Mb is the to- 
tal mass of the binary just before the onset of the CE phase 
and J is the total angular momentum of the binary at this 
point. The free parameter 7 describes the specific angular 
momentum of the ejected material (in this case the CE) in 
units of the specific angular moment um of the binary sys- 
tem. Indeed. FNelemans fc ToutI (|2005| ) found further support 
for the 7-algorithm in a few observed pre-CVs and sub-dwarf 
B p lus main sequence star binaries. 

iBeer et all l|2007l ) suggested that some systems may 
avoid the in-spiral in a CE phase altogether as a result 
of wind ejection of the envelope material from the sys- 
tem, which arises due to super-Eddington accretion on to 
the companion star. Employing a nalytical approxim ations 
for the evolution of such systems, iBeer et ahl (|2007t ) found 
that both the canonical energy budget approach and the 7- 
algorithm can describe the outcomes of the proposed mech- 
anism. 

A few observed double white dwarf binaries could also 
be explained by ap E > 1, for example, PC 1115-fl66 
ijMaxted et al.ll2002h . This has led to the suggestion that 
an energy source other than the release of orbital energy is 
being used to eject the CE from the system, such as thermal 
energy and recombination energy of ionized material within 
the giant's envelope (Han et al. 1994, Han et al. 1995, Dewi 
& Tauris 2000, W ebbink 2007). 

Alternatively, iPolitano fc Welled (|2007l ') suggested that 
rather than acE being a constant, global value for all bi- 
nary systems, it is a func tion of one of the b inary orbital 
parameters. In this spirit, iPolitano fc Welled (|2007l ) calcu- 
lated the present day population of post-common envelope 
binaries (PCEBs) and CVs if acE is a function of the sec- 
ondary mass. The investigation was motivated by the fact 
that very few, if any, CVs with brown dwarf secondaries at 
orbital periods below 77 mi ns have b e en de tected. This ap- 
pears to be in conflict with IPolitanol (|2004l ) who estimated 



from his models that such systems should make u p approx- 
imate ly 15 per cent of the current CV population. IPolitanol 
(2004") suggested that this discrepancy may be a result of the 
decreasing energy dissipation rate of orbital energy within 
the CE for decreasing secondary mass, and that below some 
cut-off mass, a CE merger would be unavoidable. 

Previous population synthesis studies into the CE phase 
of white dwarf-main sequence binaries(e.g. de Kool fc Ritter 
1993; Willems fc Kolb 2004) just considered the value of the 
envelope's binding energy, iJbiiid, due to the gravitational 
binding energy of the primary's envelope. 



, GM{r) , 
-Ebind = - / dm. 



(3) 



where M(r) is the mass contained within the radius r of 
the primary, Mc is the mass of the giant's core, and G is 
the universal gravitational constant. In population synthesis 
studies the binding energy is commonly approximated as 



-Ebind = — 



GMii\/o, 
\Ri 



(4) 



where Ri is the radius of the primary star, with a suitable 
choi ce for th e dimen sionle ss factor A. 

Ide Kooll (|1992| ) and IWiUems fc Kolbl (|2004l ) for ex- 
ample, used a constant, global value of A = 0.5. As 
iDewi fc Taurid l|2000l ) pointed out, however, this may lead 
to an overestimation of the giant envelope's binding energy 
for large radii, requiring more energy to eject it. 

In the present study we apply population synthesis tech- 
niques to calculate the present day population of PCEBs by 
considering the aforementioned theoretical descriptions of 
the CE phase. We consider different constant, global values 
of acE and A, acE as a function of secondary mass, and A 
calculated according to the internal structure of the giant 
star and the internal energy of its envelope. Also, we con- 
sider the description of the CE phase is terms of the angular 
momentum budget of the binary system. 

We then compare our model PCEB populations to 
the observed sample of white dwarf-main sequence star 
(WD-(-MS) and sub-d warf-main sequence star (sd-|-MS) 
binaries from RKCat l|Ritter fc Kolbl 120031 ), Edition 7.10 
(2008), and newly detected PCEBs from the Sloan Dig- 
ital Sky Survey (SDSS) (Rebassa-Mansergas et al. 2007; 
Rebassa-Mansergas et al. 2008). We also compare this ob- 
servational sample to our theoretical population for a range 
of initial secondary mass distributions. 

Our suite of populat ion models repre s ents a major ad- 
vance over the work by IWiUems fc Kolbl (|2004l ) who con- 
sidered only the formation rate of PCEBs, and a simplified 
estimate of their present-day populations based on birth and 
death rates. Our models also cover a more compr e hensi ve pa- 
rameter space than those bv IPolitano fc Welled l|2007l ). and 
we assess the model distributions against the location of ob- 
served systems in all three system parameters (component 
masses and period) simultaneously. 

The structure of the paper is as follows. In Section 2 
we give our computational method and models describing 
the CE phase. In Section 3 we present our results, which 
are discussed in Section 4. We conclude our investigation in 
Section 5. 
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2 COMPUTATIONAL METHOD 



Table 1 . Different models on the treatment of the CE phase. 



We employ the same method as used in iDavis" et all (|2008l ) 
to calculate the binary population, and so we present only 
a summary here. We first calculate the unweighted zero- 
age PCEB population (i.e. WD+MS systems that have just 
emerged from a CE phase) with BiSEPS (Willems & Kolb 
2002; Willems & Kolb 2004). BiSEPS er nploys the single star 
evolut ion (SSE) formulae described in iHurlev. Pols fc Tout! 
l|2000l ) an d a binary evolution scheme based on that de- 
scribed by iHurlev. Tout fc Pols! (,20021 We th en use a sec- 
ond code introduced in IWiUems et al] (|2005l ) to calculate 
the present day population of PCEBs. 

This latter step makes use of evolutionary tracks, cal- 
culated by BiSEPS, for each PCEB configuration, from the 
moment the PCEB forms to the point at which the bi- 
nary ceases to be a PCEB (for example, if the system be- 
comes semi-detached). BiSEPS therefore produces a library 
of PCEB evolutionary tracks which contain, depending on 
the population model, between 30 000 and 100 000 PCEB 
sequences. 



2.1 Initial Binary Population 

BiSEPS evolves a large number of binary systems, initially 
consisting of two zero-age main sequence stellar components. 
The stars are assumed to have a population I chemical com- 
position and the orbits are circular at all times. The initial 
primary and secondary masses are in the range 0.1 to 20 M©, 
while the initial orbital periods range from 0.1 to 100 000 d. 
There is one representative binary configuration per grid cell 
within a three-dimensional grid consisting of 60 logarithmi- 
cally spaced points in primary and secondary mass and 300 
logarithmically spaced points in orbital period. Hence we 
evolve approximately 5.4 x 10^ binaries for a maximum evo- 
lution time of 10 Gyr. For symmetry reasons only systems 
where Mi > M2 are evolved. 

The probability of a zero-age PCEB forming with a 
given white dwarf mass Mwd, secondary mass M2 and or- 
bital period Porb is determined by the probability of the bi- 
nary's initial parameters. We assume that the initial primary 
mass, Mi,i is distrib uted according to the initia l mass func- 
tion (IMF) given bv iKroupa. Tout fc Gilmord l|l993D . The 
distribution of initial secondary masses, M2,i, is obtained 
from the initial mass ratio distribution (IMRD), n(gi), where 
Qi = M2,i/Mi,i, and n{qi) has a power-law dependence on 
Qi. Alternatively, we determine M2,i from the same IMF as 
used for the initial primary mass. For brevity, we adopt 
the acronym 'IMFM2' for this latter case. Our standard 
mo del is a fiat distr i bution , i.e. n(qi) = 1. Finally, we fol - 
low llben fc Tutukovl (|1984| ) and IHurlev. Tout fc Pol j (|2002l ) 
and adopt a logarithmically flat distribution in initial orbital 
separations. 

We assume that all stars in the Galaxy are formed in 
binaries, and that the Galaxy has an age of 10 Gyr. The 
aforementioned distributions are then convolved with a con- 
stant star formation rate normalised such that one binary 
with Mi,i > 0.8 Mq is formed each year, consistent with 
observations (|Weidemannlll990l ). This gives an overall star 
formation rate of 7.6 yr~^. 
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2.2 The Common Envelope Phase 

Here we specify details of the different model descriptions 
used for the CE phase. These are also summarised in Table 
[1] In the energy budget approach the CE phase is modelled 
by equating the change in the binding energy of the pri- 
mary's envelope to the change in the orbital energy of the 
binary system (see Section 1, eqn.[l|. Thus, the flnal orbital 
separation of the binary after the CE phase, ^cE,f, is given 
by 



A 



CE,f 



Ace 



1 + 



\^ AacE»"L,i 



V M2 y 



(5) 



Here ^cE,i is the initial orbital separation at the onset of 
the CE phase, M\ and A/2 are the primary and secondary 
masses just at the start of the CE phase respectively, Mc 
is the primary's core mass, tl.i = /?L,i/^CE,i is the radius 
of the primary star in units of the orbital separation at the 
start of the CE phase and Menv is the mass of the giant's 
envelope. For our standard model we us e A = 0.5, which 
is consistent with IWiUems fc Kdbl l|2004D . This forms our 
model A. 

We consider three constant, global values of acE. For 
our reference model A we use qce = 1. We also consider 
QCE = 0.1 and 0.6. These models are denoted as CEOl and 
CE06 respectively, where the last two digits correspond to 
the value of acE^ 

Following Folitano fc Welled (|2007l ) we consider acE as 
a function of the secondary mass. We first adopt 



Qce = 



M2 



(6) 



where p is a free parameter. We consider p = 0.5, 1.0 and 
2.0. These models are denoted as PL05, PLl and PL2 re- 
spectively, where the last digits correspond to the value of 
V- 

For completeness, we also consider the second functional 
form given by 



QCE = 



A-fcut 

AI2 



M2 > Mcut 
M2 < Mcut 



(7) 



where Mcut is a cut-off mass. As lPolitano fc Weileij (|2007l ) 
we use Mcut = 0.0375, 0.075 and 0.15, which corresponds to 
xl/2, xl and x2 of the sub-stellar mass respectively. These 
models are denoted as CT0375, CT075 and CT15, where the 
digits correspond to the cut-off mass. 
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Figure 1. Values of Ag with stellar radius as calculated by the 
Eggleton code (solid line) and by the EZ code (dashed line), for 
a 6 Mq star. 

We als o obtain popula t ion m odels with values of A cal- 
culated by iDewi fc TaurisI l|2000l ). where just the gravita- 
tional binding energy of the envelope is considered (their 
Ag) and where both the gravitational binding energy and 
the thermal energy of the envelope are considered (their Ab). 
We denote these models as DTg and DTb respectively (and 
set QcE = 1 for these). 

For each of our binary configurations we calculate the 
value of Ag or Ab by applying linear i nterpolation of t he va l- 
ues for M l, Ri and Ag.b t a bulat ed bv lDewi fc TaurisI l|2000t ). 
However, iDewi fc TaurisI (|2000l ) only tabulate A values for 
primary masses M\ > 3 Mq, and for primary radii up to 
between 400 and 600 Kq . We extended this table using full 
stellar models calculated by the Eggleton code (provided by 
van der Sluys, priv. comm.), up to the point where the dy- 
namical timescale of the star becomes too small, and the 
evolutionary calculations stop. 

In particular, for the 1 Mq stellar model, we obtain 
values of Ag and Ab for radii up to approximately 500 Rq. 
For stellar masses between 2 and 4 Mq we obtain A values up 
to approximately 1000 R©, while for stellar masses between 
5 and 8 MeQ up to between approximately 400 and 600 R© . 

However, we find that PCEB progenitors with primary 
masses between 2 and 8 Mq can fill their Roche lobes with 
radii up to about 1200 Kq. To estimate likely A values for 
these large radii we used the EZ stellar evolution code, which 
is a r e- written, stri pped down adaptation of the Eggleton 
code (|Paxton ' 2004). Indeed, the code was designed to be 
more stable and faster, though at the cost of some detail. In 
the mass range of approximately 5 to 6 Mq this does return 
Ag values for large radijfl (see Fig. [1] for the 6 Mq star). 
However, for masses below 5 Mq or above 6 Mq, the EZ code 
stops the evolution at radii similar to the Eggleton code. 
For these stellar masses, we simply take the last calculated 
values of Ag for any _Ri > 600 Rq in our population models. 

If there is sufficient thermal energy within the giant's 
envelope such that its total binding energy becomes positive 
(Ab < 0) the formalism described by equation ((SJ breaks 
down. We model this instead as an instantaneous ejection of 
the giant's envelope via a wind, which takes away the specific 

^ The upper limit of 8 Mq is the largest primary progenitor which 
will form a white dwarf. 

^ For these masses, Ab < at radii Fti > 500. 



orbital angular momentum of the giant. What remains is a 
wide WD-I-MS binary. We do not consider these any further 
in this investigation. 

2.3 CE Phase in terms of Angular Momentum 

WefoUow iNelemans et all (|2000l ) and iNelemans fc Tout! 
l|2005l ) by considering models where the CE phase is de- 
scribed in terms of the change of the binary's angular mo- 
mentum, AJ, given by equation in Section 1. If J and Jf 
are the total angular momentum of the binary immediately 
before and after the CE phase respectively, then equation 
^ can be re-written as 

At the start of the CE phase, we not only consider the or- 
bital angular momentum of the binary, Joib.i, but also the 
spin angular momentum of the Roche lobe-filling giant star, 
Jspin.i, such that J = Japin.i + Jorb.i- If the giaut is rotating 
synchronously with the orbit with a spin angular speed of 
f2orb, then 

Jspin,! = fc^Mlii^^i^orb, (9) 

where is the ratio of the gyration radius to the radius 
of the giant star. W e use the values of fc^ tabulated by 
iDewi fc TaurisI (200&). Our procedure for calculating for 
each giant star configuration is the same as that used to cal- 
culate Ag and Ab. The initial orbital angular momentum of 
the binary, on the other hand is given by 

Similarly the final orbital angular momentum of the binary 
after the CE phase, Jorb,f , is given by 

Note that we do not consider the spin angular momentum of 
either of the stellar components after the CE phase as these 
are negligible compared to the orbital angular momentum 
of the binary. Thus, Jf = Jorb,f- From equations (|9} to (|11|) 
we can therefore solve for AcE.f in equation (|8]). 

For the large majo rity of the observed systems that 
iNelemans fc Tout! (|2005l ) considered, they found that their 
CE phase could be reproduced with 7 « 1.5. We therefore 
consider 7 = 1.5 in equation ([Sjl. This is our model nl5, 
which is summarised in Table [1] 

2.4 Magnetic Braking 

After the CE phase, the subsequent evolution of the PCEB 
will be driven by angular momentum losses from the or- 
bit. We follow the canonical disrupted magnetic braking hy- 
pothesis, and assume that for a fully convective secondary 
star with M2 < Mconv.MS = 0.35 Mq (where Mconv,MS is the 
maximum mass of a fully convective, isolated, zero-age main 
sequence star) , gravitational radiation is the only sink of an- 
gular momentum. For secondaries with M2 > Mconv,MS on 
the other hand, we assume that the evolution of the PCEB 
is driven by a combination of gravitational radiation and 
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magnetic braking. For this investigation we consider the an- 
gular mo mentum loss rate due to ma gnetic braking, Jmb, as 
given bv iHurlev! Tout fc Polj (|2002[ ) 



(12) 



where R2 is the secondary's radius, Mconv is the mass of the 
secondary's convective envelope and is its spi n frequenc y. 
We n ormalise equation (|12p as described in i Davis et al . 
l|2008l ) by applying a factor = 0.19 which gives the angu- 
lar momentum loss rate appropriate for a period gap of the 
observed width of one hour. Magnetic braking also becomes 
ineffective for PCEBs with M2 > 1.25 M©, which have very 
thin or no convective envelopes. 



3 RESULTS 

We begin by comparing our theoretical population of PCEBs 
from model A, n{qi) — 1, with observed WD-I-MS and 
sd-l-MS syst ems, obtained from Edition 7.10 (2008) of 
RKCat (Rit ter fc Kolbl l2003t ). This sample is also supple- 
mented with newly discovered PCEBs from SDSS (Rebassa- 
Mansergas et al. 2007; Rebassa-Mansergas et al. 2008). Ob- 
served WD-I-MS systems are tabulated in Table (2] while 
observed sd-|-MS systems are tabulated in Table [S] We then 
discuss the present day number and local space densities of 
PCEBs, and discuss each of the CE models in turn. 



3.1 White dwarf- main sequence systems 

Figure [2] shows the theoretical present-day PCEB popula- 
tions on the M2 — log Porb plane for model A, assuming an 
IMRD of n{qi) = 1, along with the observed WD-f MS sys- 
tems. 

Both the theoretical distributions and observed systems 
are divided into panels according to their white dwarf masses 
as indicated in each of the nine panels, labelled (a) to (i). The 
grey-scale bar at the top of the plot indicates the number of 
systems per bin area. We take into account the uncertainty 
of the white dwarf masses for each observed system by cal- 
culating the weighting of the system in each panel that over- 
laps with the white dwarf uncertainty interval. The panels 
in which these systems lie, and the associated weightings, 
are also shown in Table [2] The weightings are calculated 
from a Gaussian distribution with a mean white dwarf mass 
(Mwd) and a standard deviation a, which corresponds to 
the measured white dwarf mass and its uncertainty respec- 
tively. The colour bar on the right hand side of the plot 
indicates the weightings. 

There is an acceptable overall agreement between the 
observed systems and the theoretical distributions, in the 
sense that all but two observed systems are in areas popu- 
lated by the standard model. However, the model distribu- 
tions extend to areas at orbital periods longer than about 1 
d and donor masses larger than about 0.5 M0 with very few, 
if any, observed systems. The two outliers that cannot be ac- 
counted for by our standard model are V651 Mon (Mendez 
& Niemela 1981; Mendez et al. 1985; Kato, Nogami & Baba 
2001) and IK Peg (Venues, Christian & Thorstensen 1998). 
The location of these systems are indicated in Figure [3] by 
the arrows (these systems lie outside of the range displayed 



in Figure [2] and so are not shown there) . We shall discuss 
these systems in turn to deduce whether they are in fact 
PCEBs. 

IK Peg contains a 1.2 Mq white dwarf, and a 1.7 Mq 
secondary star with spectral type A8 (Landsman, Simon 
& Bergeron 1993; S malley et al. 1996; V enues, Christian 
& Thorstensen 1998). ISmaUev et al] l| 19961 ) found that the 
secondary star had an overabundance of iron, barium and 
strontium, which m ay be accounted for if IK Peg under- 
went mass transfer. ISmallev et al.l l|l996l) suggest that this 
was in the form of a CE phase. As the secondary plunged 
into the envelope of the giant it would be contaminated by 
s-processed material. As WD-I-MS systems can also form 
through a thermally unstable Roche lobe overflow (RLOF) 
phase (see Willems & Kolb 2004, formation channels 1 to 3), 
we calculated the present day population of WD-I-MS sys- 
tems with 1.1 < A'/wd/Mq < 1.44 that form through such 
a case B thermal-timescale mass transfer (TTMT) phase 
(channel 2 in Willems & Kolb 2004), which is shown in the 
right panel of Figure |4] The diamond indicates the location 
of IK Peg. As the location of IK Peg cannot be accounted for 
by our theoretical population of case B TTMT systems, it is 
hence a likely PCEB candidate. We explore this possibility 
further in Section 3.5 where we consider more population 
synthesis models of the CE phase. 

V651 Mon contains a 0.4 Mq hot white dwarf primary 
with an A-type secondary star. This system lies within a 
planetary nebula, which cou l d be th e remnant of a CE phase. 
However, Ide Kool fc RitteJ l|l993l) suggest that the plane- 
tary nebula may have formed when the compact primary 
underwent a shell flash, and instead the binary formed via 
a thermally unstable RLOF phase. As for IK Peg, we cal- 
culate the present day population of WD-fMS systems with 
A^wd/M© < 0.4 that formed through a thermally unsta- 
ble RLOF phase, shown in the left panel of Figure U) This is 
compared with the location of V651 Mon shown as indicated 
by the arrow. From Figure |4] it appears hence feasible that 
V651 Mon formed from a case B RLOF phase. 

To obtain the possible progenitor of V651 Mon, we first 
find the calculated WD-I-MS configuration which formed 
through the case B TTMT channel which is situated closest 
to the location of V651 Mon in {Mwo, M2, Porh) space. If 
5 Mwd, SM2 and SPoih is the difference between the observed 
white dwarf mass A/wD.obs, secondary mass Afa.obs and or- 
bital period Porb.obs of V651 Mon and a calculated WD+MS 
binary configuration, then the normalised distance. A, in 
(AfwD, A^ajPorb) space is given by 



A = 





<5AfwD 


2 

4- 


5M2 


2 
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SPoih 


1 




AfwD,obs 


Af2,obs 


p3rb,obs 





1/2 



(13) 



Thus we require the WD-I-MS binary configuration that 
gives the smallest value of A. We find that A is minimised 
for A^wD = 0.33 Mq, M2 = 1.76 M© and Poa = 15.92 d. 
This system has a ZAMS progenitor with Mi,i = 2.47 Mq, 
M2,i = 0.98 Mq and Po,b,i = 2.60 d. 

While the secondary mass of this WD-I-MS configura- 
tion is within the measured uncertainty of the secondary 
mass in V651 Mon, its white dwarf mass is not, and in 
fact underestima t es th e wh ite dwarf mass of V651 Mon. 
iTout fc EggletonI l| 19881 ) and lTout fc EggletonI (|l988bl ) sug- 
gest that V651 Mon may have formed with enhanced mass 
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Table 2. The orbital periods, white dwarf masses (Mvvd), and secondary masses, M2, of the observed sample of WD+MS systems. 
Also shown are the weightings of each system in the associated white dwarf mass range, corresponding to panels (a) to (i) (see Fig. 
0. 



Name 


Alt. name 




/Mq 


M2 

/Mq 


Panel 


Weighting (%) 


Ref. 


V651 Mon 


AGK3-0 965 


15.991 


0.40+0.05 


1.80+0.3 


a,b 


50.0,48.0 


1,2 


SDSS J1529+0020 




0.165 


0.40+0.04 


0.25+0.12 


a,b 


50.0,49.0 


3 


CC Get 


PC 0308+096 


0.287 


0.40+0.11 


0.18+0.05 


a,b,c 


50.0,32.0,1.00 


4 


LM Com 


CBS 60 


0.259 


0.35+0.03 


0.17+0.02 


a 


95.0 


5 


HR Cam 


WD 0710+741 


0.103 


0.41+0.01 


0.096+0.004 


a,b 


16.0,84.0 


6,7 


0137-3457 




0.080 


0.39+0.035 


0.053+0.006 


a,b 


61.0,39.0 


8 


SDSS J2339-0020 




0.656 


0.8+0.4 


0.32+0.09 


a,b,c,d 


14.0,7.0,4.0,4.0 


3 












e,f,g,h,i 


5.0,5.0,10.0,27.0,17.0 




SDSS J1724+5620 




0.333 


0.42+0.01 


0.25 - 0.38 


b 


98.0 


3 


RR Cae 


LET 349 


0.303 


0.44+0.023 


0.18+0.01 


b 


95.0 


9 


MS Peg 


CD 245 


0.174 


0.49+0.04 


0.19+0.02 


b,c 


59.0,33.0 


10 


HZ 9 




0.564 


0.51+0.1 


0.28+0.04 


b,c,d,e 


32.0,10.0,16.0,10.0 


11,12,13 


GK Vir 


PC 1413+015 


0.344 


0.51+0.04 


0.10 


b,c 


40.0,44.0 


14,15 


LTT 560 




0.148 


0.52+0.12 


0.19+0.05 


b,c,d,e 


28.0,16.0,15.0,11.0 


16 


DE CVn 


J1326+4532 


0.364 


0.54+0.04 


0.41+0.06 


c,d 


44.0,33.0 


17 


2237+8154 




0.124 


0.57+0.1 


0.30+0.1 


b,c,d,e,f 


20.0,18.0,20.0,17.0,12.0 


18 


SDSS J1151-0007 




0.142 


0.6+0.1 


0.19+0.08 


b,c,d,e,f 


14.0,15.0,19.0,19.0,15.0 


3 


UZ Sex 


1026+0014 


0.597 


0.68+0.23 


0.22+0.05 


b,c,d,e. 


11.0,6.90,7.80,8.40, 


4 












f,g,h 


8.7,16.0,27.0 




NN Ser 


PC 1550+131 


0.130 


0.54+0.04 


0.150+0.008 


c,d 


44.0,33.0 


19 


FS Get 


Eeige 24 


4.232 


0.57+0.03 


0.39+0.02 


c,d 


24.0,59.0 


20,21 


J2013+4002 




0.706 


0.56+0.03 


0.23+0.01 


c,d 


35.0,54.0 


21,22 


J2130+4710 




0.521 


0.554+0.017 


0.555+0.023 


c,d 


41.0,59.0 


23 


1042-6902 




0.337 


0.56+0.05 


0.14+0.01 


c,d,e 


31.0,37.0,18.0 


21 


SDSS J0314-0111 




0.263 


0.65+0.1 


0.32+0.09 


c,d,e,f,g 


9.0,15.0,19.0,19.0,24.0 


3 


IN CMa 


J0720-3146 


1.262 


0.58+0.03 


0.43+0.03 


d,e 


59.0,24.0 


24 


J1016-0520AB 




0.789 


0.61+0.06 


0.15+0.02 


d,e,f 


28.0,31.0,19.0 


24 


EG UMa 


Case 1 


0.668 


0.63+0.05 


0.36+0.04 


d,e,f 


22.0,38.0,26.0 


25 


2009-1-6216 




0.741 


0.62+0.02 


0.189+0.004 


e 


77.0 


26 


BE UMa 


PC 1155+492 


2.291 


0.70+0.07 


0.36+0.07 


e,f,g 


16.0,26.0,42.0 


27,28 


1857-f5144 




0.266 


0.80+0.2 


0.23+0.08 


e,f,g,h 


6.80,8.20,19.0,43.0 


29 


SDSS J0246+0041 




0.726 


0.9+0.2 


0.38+0.07 


f,g,h 


5.0, 15.0,53.0 


3 


SDSS J0052-0053 




0.114 


1.2+0.4 


0.32+0.09 


g,h,i 


5.0,24.0,32.0 


3 


BPM 71214 




0.202 


0.77+0.06 


0.540: 


g,h 


57.0,31.0 


30 


QS Vir 


1347-1258 


0.151 


0.78+0.04 


0.43+0.04 


g,h 


67.0,31.0 


31 


V471 Tau 


BD+16 516 


0.521 


0.84+0.05 


0.93+0.07 


g,h 


21.0,79.0 


32 


IK Peg 


BD+18 4794 


21.72 


1.19+0.05 


1.7 


i 


96.0 


33 



Re ferences: (l)|M endcz & Nicm elal lll981^. (2)lMendez et aU (Il985l). (3 )|Rcbassa-Manser gas et al.l ||2008|) . (4')ISaffer et al.l ( 1993 ). 
('5') IShimanskv e? al., C200y . (•6l iMarsh fc Ducld 19961 ). (7llMaxted et al.l lll99a')._ ('8) ,Ma.xted et al.l ll2006l\. ('9')lBragaglia et al.l 1 199E 
(lOllSchmidt et al.l lll995l) 
lll993ll. fl5')lGreen et~aLl 



Stauffej ||19'87F 



12) | Lanning . 

■ _. ■ II1978I'). fl6')lTaDpert et al.l (l2007l 'l". QTlh 

ri9llCatalan et al.l ||1991k (2 0')IVerines fc ThorstensenI (Il994h, f2 1') lKawka et al.l 1I2OO8I'). f22')lGood et al.l |20^7l23)]M^te^^] 



Pesch 



1200 



411. ('24)IVennes et al.l 



(28) 
(32) 



Fergusor^_et_^ 
O'Brien et al.l 



.^^^^ .,199 9!). (25VBlea ch et al.l (|2000D . (26)'Morales-RuedaetaL 
l[l999h . (291 Amigwero iwit et al.l ||2007^, C30') .Kawka ct al] 1I2OO2I) 
200lh . ('33') lLandsman. Simon fc BergeronI l|l993h 



lll98j). (lillGuinan fc SionI lll984h. (14llFulbri'ght' et al.l 
den Besselaaj ll2007h. (18)lGansicke et al.l 
od et al.] | |2005^, ('23')lMax 
i|2005l ) . r27)|Wood et al 



f31') IO'Donoghue et al 



1995 ). 



20031'! . 



loss from the primary star before it filled its Roche lobe. 
Thus, a more massive progenitor primary than the one found 
in our best-fit WD+MS configuration may have formed the 
white dwarf. If enough mass is lost from the primary due 
to enhanced wind losses before it fills its Roche lobe to suf- 
ficiently lower the mass ratio, M-1/M2, then the resulting 
mass transfer will be dynamically stable. 

3.2 Sub-dwarf + main sequence star binaries 

One of the main formation chann els of sd+MS binaries is 
via the CE phase l|Han et al.ll2002l ): if the primary star fills 
its Roche lobe near the tip of the first giant branch, and the 



CE is ejected from the system, then a short-period binary 
with an sd star is formed once the helium core is ignited. 
iHan et al.l (|2002l ) found that the mass of the sdB star ranges 
from approximately 0.32 to 0.48 Mq. The typical mass for 
a sdB star is approximately 0.46 M©. This appears to be 
broadly consistent with the observed masses of sd stars (see 
Table |3}- The mass of the sub-dwarf in V1379 A ql is some- 
what smaller than the calculated lower limit by iHan et al.l 
(l2002t ). while the sub-dwarfs in UU Sge and KV Kel are 
more massive than the calculated upper limit in sub-dwarf 
mass. 

The formation of sub-dwarf stars is affected by the wind 
loss of the giant's envelope, the metallicity of the primary, 
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Table 3. The orbital periods, sub-dwarf mass (Msd)i secondary mass, M2, and spectral type, SpT, of the sub-dwarf for the 
observed sample of sd+MS systems. Also shown are the weightings of each system in the corresponding panels (a) to (f) (see Fig. 
[5)l. except for those systems where no uncertainty in the sub-dwarf mass is available. A colon denotes an uncertainty flag. 



Name 


Alt. name 


Porb/ d 


AfsD 


SpT 1 


M2 


Panel 


Weighting (%) 


Ref. 








/Mq 




/Mq 








V1379 Aql 


HD 185510 


20.662 


0.304±0.015 


sdB 


2.27±0.13 


a 


100.0 


1 


FF Aqr 


BD -3 5357 


9.208 


0.35±0.06 : 


sdOB 


1.40±0.4 : 


a,b 


80.0,20.0 


2,3 


2333-1-3927 




0.172 


0.38±0.09 


sdB 


0.28±0.04 


a,b 


59.0,32.0 


4 


AA Dor 


LB 3459 


0.262 


0.330±0.006 


sdO 


0.066±0.001 


a 


100.0 


5 


HW Vir 


BD -7 3477 


0.117 


0.48±0.09 


sdB 


0.14±0.02 


a,b,c,d 


19.0,40.0,19.0,13.0 


6 


2231-1-2441 




0.111 


0.47 : 


sdB 


0.075 


b 




7 


NY Vir 


PG 1336-018 


0.101 


0.466±0.006 


sdB 


0.122±0.001 


b 


100.0 


8 


0705-1-6700 




0.096 


0.483 


sdB 


0.134 


b 




9,10 


BUL-SC 16 335 




0.125 


0.5 : 


sdB 


0.16 


b 




11 


J2020-I-0437 




0.110 


0.46 


sdB 


0.21 


b 




12 


XY Sex 


1017-0838 


0.073 


0.50 


sdB 


0.078±0.006 


b 




13 


V664 Cas 


HFG 1 


0.582 


0.57±0.03 


sdB 


1.09±0.07 


c,d 


24.0,59.0 


14,15,16 


UU Sge 




0.465 


0.63±0.06 


sdO 


0.29±0.04 


d,e,f 


22.0,32.0,25.0 


17 


KV Vel 


LSS 2018 


0.357 


0.63±0.03 


sdO 


0.23±0.01 


d,e,f 


15.0,59.0,24.0 


18 



Refere nces : (DlJefferv fc SimorJ lll997t). f2')IVaccaro fc WilsonI ||2003|). (3^1261 et al.l l ll977l). (4')lHeber et al.) (l2o3ki5}|Rauch 
1 20001') . C6')IWood fc Saffcri 1*1999''). ('7V0ste nsen et al.lll2007^ . (S llKilkennv et al.l ll l998"). (9VDrcchscl ct al.' ('2001'). aO')'Ncmcth et al. 
1 2005 ). (irj TPolubck et al.,, 12007 ). ('12) Wils et al.l (|2007|) . (1 3')lMaxted et al.rl20 02). (14) Miller ct al.. 11976) . (15) ,Pollacco fc Bell 
I I993I ). a6'i lBell et al l il994 ). fl7) lHilditch et al.l (119961) ? (18) IShimanskii et al.lT2004i ) 



and the degree of convective overshooting of the primary. A 
detailed consideration of sdB binaries is beyond the scope 
of this study. Instead, we simply plot the observed sample 
of sd-|-MS binaries over our calculated distributions of the 
■WD-I-MS PCEBs (see Figure [5] for model A and n{qi) = 1, 
in the same style as Figure [2}, thereby assuming that there 
is little orbital evolution of these systems by the time the 
sdB stars become white dwarfs. This is justified as follows: 

The evolution of sd-l-MS binaries will be driven by wind 
losses from the subdwarf pri mary with a mass loss rate of a p- 
proximately 10"" Mg yr'^ (|Vinkl (|2004h : IUnglaubl (|2008l )'). 
as well as systemic angular momentum losses. For the ob- 
served sample of sd-|-MS systems in Table |31 M2 < 0.35 Mq 
in all cases, and so gravitational wave radiation will there- 
fore be the only sink of orbital angular momentum. During 
the lifetime of an sdB (i.e. t he core helium burning tim e) of 
approximately few x 10* yr (jHurlev. Pols fc Toutll200oh . we 
estimate a relative change in the binary's orbital period of 
less than about 7 per cent. 

Figure [5] shows that, as with the observed 'WD-I-MS 
systems, there are observed sd-|-MS systems that cannot be 
explained by our standard model. These are FF Aqr and 
V1379 Aql. Their locations are indicated in the left panel of 
Figure O by the arrows. 'We discuss these systems individu- 
ally to determine whether they are in fact PCEBs. 

VI 379 Aql This system contains a K-giant secondary 
and an sdB primary, which may be a nascent helium 
white dwarf |jefferv. Simon fc EvansI [l992l ) . Furthermore, 
IJefferv fc SimonI (Il997l ) argue tharV1379 Aql formed from 
a thermally unstable RLOF phase. As a consequence of the 
synchronous rotation of the secondary, the system exhibits 
RS CVn chromospheric activity due to the K-giant 's en- 
hanced magnetic activity. This is likely to occur if the system 
had a previous episode of RLOF, rather than if the system 
has emerged from a CE phase. 

iTout fc EggletonI l|l988bl ) suggest that the progenitor 
binary of 'V1379 Aql had Mi,i = 2.0 Mq and M2.i = 1.6 



Mq. The secondary star would have completed one half 
of its main sequence lifetime by the time the primary 
evolved to become a giant aft er approximately 6 x 10* yr 
(iJefferv. Simon fc Evan'3 Il992l '). If a CE phase did occur, 
then the secondary would still have been on the main se- 
quence, and hence possess an insubstantial convective enve- 
lope. Consequently, magnetic braking would be ineffective 
in shrinking the orbital separation sufficiently for tidal in- 
teractions to enforce synchronous rotation of the secondary 
star. 

iNelemans fc ToutI (|2005l ). in reconstructing a possible 
CE phase for this system, could not find a solution for acE, 
and cited this as evidence for their alternative CE descrip- 
tion in terms of the angular momentum balance. We suggest 
that no solution was found because of the possibility that 
V1379 Aql formed from a thermally unstable RLOF phase. 

The left panel of Figure |4] compares the theoreti- 
cal 'WD-I-MS systems with A/wd/Mq < 0.4 that formed 
through a thermally unstable RLOF phase with the loca- 
tion of V1379 Aql shown by the arrow. It appears hence 
feasible that 'V1379 Aql formed through this evolutionary 
channel. As with 'V651 Mon, we apply equation (|13|l to find 
the 'WD+MS configuration which lies the closest to the lo- 
cation of 'V1379 Aql in (AfwD, M2, Porb). Doing this we find 
a 'WD-fMS configuration with Mwd = 0.23 Mq, M2 = 2.31 
M0 and Porb = 20.58 d. This system forms from a ZAMS 
progenitor with A'/i,i = 1.79 Mq, M2,i = 0.77 Mq and 
Porb.i = 1.41 d. 

The mass of the white dwarf in our best-fit 'WD-I-MS 
configuration onc e again underestimates th at of Y1379 Aql. 
As for V651 Mon lTout fc EggletonI (|l988bh suggest that the 
primary star in V1379 Aql also underwent enhanced mass 
loss before it filled its Roche lobe. Thus, the mass ratio, 
M1/M2, was decreased sufficiently so that a dynamically 
unstable RLOF phase was avoided. The actual progenitor 
primary of V1379 Aql may therefore have been more massive 
than what we have predicted in our calculations. Note also 
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Figure 2. Calculated present-day PCEB populations, for model A and for n{q{) = 1. The distributions over A/2 and logPo^b are shown 
for nine different white dwarf mass intervals as indicated in each of the nine panels, labelled (a) to (i). The grey-scale bar at the top of 
the plot in dicates the number o f systems per bin area. Also shown are the known WD-I-MS systems obtained from Edition 7.10 (2008) 
of RKCat l|Ritter fc Kolbll2003h . We take into account the uncertainties of the white dwarf mass for each system by calculating the 
weighting of each system within each panel. This is repeated for those systems whose error bars in white dwarf mass overlaps more 
than one panel. The weighting is indicated by the colour bar on the right hand side of the plot. The dotted line in panel (b) indicates 
that the secondary mass in SDSS J1724-I-5620 lies between 0.25 Mq and 0.38 Mq. 



that in our calculations, we do not treat the sdB phase of 
the white dwarf explicitly. 

FF Aqr also dis plays RS CVn characteristics 
(IVaccaro fc WilsonI l2003l i and hence we can apply the 
same argument used for V1379 Aql to FF Aqr; the 
secondary has a mass of approximately 1.4 Mq and so 
will have an insubsta ntial convective envelope. Indeed, 
Ide Kool fc Ritteij (liggsl ) argue that FF Aqr formed from a 
thermally u nstable RLOF phase, p erhaps with enhanced 
wind losses (|Tout fc Eggletonlll988bl ). 

The left panel of Figure |4j compares the theoretical 
WD-fMS binary population with Mwd /Mq < 0.4 , which 
formed through a channel 1 IWillems fc Kolbl (|2004l ') RLOF 
phase, with the location of FF Aqr indicated by the arrow. 
This supports the possibility that FF Aqr formed from a 
TTMT RLOF phase. While our best fit WD-^MS configu- 
ration has A/wD = 0.30 M0, which is consistent with the 
observed value within the measured uncertainty, our best- 



fit configuration slightly underestimates the secondary mass 
with Ah = 1.22 Mq. 



3.3 The shape of the distributions 

We will now discuss common features of the M2 — logioPorb 
PCEB distributions, with the middle panel of Figure [3] as a 
typical example. 

Figure [6] shows the upper boundary in the population 
of progenitor systems at the onset of the CE phase which 
will form PCEBs with 0.4 < Mwd/Mq < 0.5 (long dashed 
line). The short dashed line shows the corresponding upper 



^ We note that the shapes of the di stributions for model nl5 are 
very different. We refer the reader to lNelemans fc ToutI ll2005h for 
a detailed discussion. 
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Figure 3. The theoretical present day PCEB population for model A, = 1, for a wider range in M2 and i^orb 

than Fig. [2 for 

selected white dwarf mass intervals. Each panel labelled (a) to (c) represents the indicated range in white dwarf or sub-dwarf mass. 
Diamonds: WD+MS systems; triangles: sd+MS systems. The location of the individual systems V651 Mon, V1379 Aql and FF Aqr 
are indicated by the arrows. The grey-scale on the top of the plot indicates the number of systems per bin area, while the bar on the 
right hand side indicates the weighting of the systems in their corresponding panels. 
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Figure 4. Left panel: The theoretical population of WD-I-MS systems with Mwd/Mq < 0.4 that formed through a thermally unstable 
RLOF phase, compared with the location of FF Aqr, V651 Mon, and V1379 Aql, as indicated by the arrows. Right panel: The theoretical 
population of WD+MS systems with 1.1 < Mwo/Mq < 1.44 that formed through a thermally unstable RLOF phase with a naked 
helium star remnant. This is compared with the location of IK Peg indicated by the arrow. A triangular symbol indicates that the 
system is a sd-|-MS binary, while a diamond indicates the system is a WD+MS binary. 



boundary in the resulting population of PCEBs (this repre- 
sents the upper boundary of the distribution shown in the 
middle panel of Figure [HJ, while the thick, black line shows 
the orbital period at which the secondary star will undergo 
Roche lobe overflow (RLOF), and hence the PCEB will be- 
come semi-detached. 

The long dashed boundary in Figure[S]arises because for 
each mass M2 there is an upper limit for the orbital period of 
the progenitor binary that will still form a white dwarf with 
a mass within the considered mass interval. Longer-period 
systems will either form a more massive white dwarf, or 
remain detached and not undergo a CE phase. More specif- 



ically, the upper boundary corresponds to the most massive 
white dwarfs in the panel. In the example shown this is 0.47 
M0 (not 0.50 M0) because there is a gap in the white dwarf 
mass spectrum between the low-mass helium white dwarfs 
and t he more massive carbon-ox ygen white dwarfs. The SSE 
code (|Hurlev. Pols fc Toutll2000l ) places this gap in the range 
between 0.47 < Mwd/M© < 0.6. 

This long dashed upper pre-CE boundary maps on to 
the short dashed PCEB boundary, also shown in Figure |6l 
The masses of the progenitor primaries at the start of the CE 
phase are indicated along the long-dashed boundary, while 
the white dwarf masses are shown along the short-dashed 
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Figure 5. Same as Figure [2] but this time showing the observed sample of sub-dwarf+MS binaries. The black points represent those 
systems for which we have not calculated the weightings, as measured uncertainties for the sub-dwarf masses are not available. 



boundary of the PCEB population, at the corresponding 
point along the boundary. 

For a PCEB to form, the ZAMS progenitor primary 
must fill its Roche lobe and commence the CE phase be- 
fore the secondary itself evolves off the main sequence. This 
point is of little consequence for M2,i ^ 1 M© as the main 
sequence lifetime of such secondaries will be rMS,2 ^ 10 Gyr, 
which is on the order of the Galactic lifetime (and the max- 
imum evolution time we consider). Thus, in such cases, the 
progenitor primary must ascend the giant branch and fill its 
Roche lobe within the lifetime of the Galaxy. 

We find that the least massive ZAMS progenitor which 
will subsequently form a 0.47 M© white dwarf within the life- 
time of the Galaxy, is 1.1 Mq. Hence, for the reasons stated 
above, M2.i < 1.1 M0, if the evolution of the secondary off 
the main sequence is to remain unimportant. As indicated 
along the long-dashed boundary in Figure |S] for M2 < 1.1 
M0(to the left of the dashed vertical line), the progenitor 
primaries, as a result of wind losses, commence the CE phase 
with masses of approximately 0.9 Mq. Furthermore, as the 
long-dashed boundary for M2 < 1.1 Mq is approximately 
flat, the progenitor primaries along this boundary fill their 
Roche lobes at the same orbital period. 

To see how this determines the shape of the resulting 
short-dashed PCEB boundary for M2 < 1.1 Mq (labelled 
'A' in Figure [6]) we consider equation ([Sj, which shows that 
for constant Mc (and hence the white dwarf mass). Mi and 
rL,i (at the start of the CE phase), Ace,i decreases for de- 
creasing M2. Hence the orbital period of the short-dashed 



boundary to the left of the dashed line in Figure [6] decreases 
for decreasing A'h- 

On the other hand, for M2,i > 1.1 Mq (to the right of 
the dashed line in Figure |6} , the main sequence lifetime of 
such stars is less than 10 Gyr, and so the progenitor pri- 
mary must increasingly compete against the secondary to 
evolve off the main sequence flrst. As a result, the ZAMS 
progenitor primary star needs to be increasingly more mas- 
sive for increasing secondary mass. This is reflected in the 
masses of the primaries at the start of the CE phase, as 
indicated along the long-dashed boundary in Figure [§] for 
M2 > 1.1 Mq, which increase from approximately 0.9 Mq 
where M2 « 1.1 Mq, to 1.8 Mq where M2 « 1.8 Mq. 

For a primary star ascending the giant branch with lu- 
minosity Li, its radius on the giant branch, Rqb, can be 



modelled as 
1.1 



Rgb 



(Lr + 0.383Lr') 



(14) 



(|Hurlev. Pols fc Tout|[2000l '). where Lioc Mc^. Thus the ra- 
dius of the star on the giant branch is a function of its core 
mass and (weakly) of its initial ZAMS mass. For a given core 
mass, in this case of Mc = 0.47 Mq, the corresponding value 
of Rgb will decrease for increasing mass of the ZAMS pri- 
mary. For Roche lobe filling stars we have Porb oc R^^^ /M^^^ 
so the orbital period at which the CE phase will occur will 
decrease with increasing M2. This is shown as the slope in 
the long-dashed boundary in Figure [S] for M2 > 1.1 Mq. 

The shape of the short-dashed boundary to the right of 
the dashed line in Figure |6] (labelled 'B') is a consequence of 
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Figure 6. Critical boundaries in the \o^^Q(Pa-c\^) — M2 plane. The 
long dashed line shows the upper boundary in the population of 
progenitor binaries at the start of the CE phase which will form 
PCEBs with 0.4 < Mwd/Mq < 0.5. The values along this hne 
denote the mass of the Roche lobe-filling progenitor primary at 
the corresponding locations along the boundary. The short dashed 
line is the upper boundary in the resulting population of PCEBs 
with 0.4 < MwD < 0.5. The values along this line denote the 
white dwarf mass The label 'A' corresponds to the portion of the 
boundary at M2 < 1.1 Mq (short-dashed line), while 'B' denotes 
the portion of the boundary at M2>1.1 Mq. The thick black 
line is the Roche lobe-overflow (RLOF) limit, i.e. the point at 
which the PCEB will become semi-detached, and hence defines 
the lower limit in the distribution. Finally, the dot-dashed line 
corresponds to the cut-off in this PCEB population at M2 ^ 1.8 
Mq. 

the shape of the corresponding portion of the long-dashed 
pre-CE boundary. Furthermore, for increasing M2, while Mc 
remains constant, A/i,i does increase, and so the mass of the 
primary giant's envelope will correspondingly increase. The 
result of this is that AcE,f will decrease for increasing A/2, 
as equation ([5| shows. This explains why the short-dashed 
PCEB boundary to the right of the dashed lin Figure [6] 
decreasing for increasing M2. 

Note that the cut-off in the population, as indicated by 
the vertical dot-dashed line (at A'/2 ~ 1-8 M0) in Figure|6] is 
a consequence of the fact that the most massive progenitor 
primary which will form a white dwarf in the range 0.4 < 
A/wd/M© < 0.5 is 1.8 Mq. For PCEB distributions with 
larger white dwarf masses, which will be formed from more 
massive primary progenitors, this cut-off will shift towards 
larger values of Af2. 

3.4 The present day population and space 
densities of PCEBs 

Table |4] summarises the formation rate, present day popula- 
tion and the local space density, q, of PCEBs for each model 
describing the treatment of the CE phase, and for each pop- 
ulation model we considered. We calculate the local space 



density of PCEBs by dividing the present day population by 
the Galactic volume of 5 x 10^^ pc"^. 

For a given initial secondary mass distribution, we see 
that the present day population of PCEBs increases with 
increasing values of acE. By increasing qce, less orbital en- 
ergy is required to eject the CE and therefore the binary 
system will undergo less spiral in. Hence, fewer systems will 
merge during the CE phase. For a given IMRD, we find very 
little difference in the calculated formation rates, present day 
numbers and local space densities of PCEBs between mod- 
els A and DTg. For model DTb, on the other hand, where 
we are now considering the internal energy of the giant's 
envelope, we obtain a modest increase in the present day 
population of PCEBs compared to model A. 

The present day population of PCEBs decreases slightly 
with increasing va lues of p in equat ion (|6]). As shown in 
Figure 1 of Polita no fc Weileij (I2007D the average value of 
acE for a given range in Af2 decreases with increasing p. 
By increasing p more orbital energy is required to expel the 
envelope from the system, hence a larger spiral-in of the 
binary. This will in turn lead to more mergers. 

On the other hand, the present day number (local space 
density) of PCEBs decreases from model CT0375 to model 
CT15. As we increase the value of A/cut the typical value 
of Qce for a given range of Mi will decrease. Hence, fewer 
systems will survive the CE phase. 

We notice that the model nl5 which describes the CE 
phase in terms of the binary's angular momentum gives the 
largest present day population of PCEBs. For n(ci\) = 1, 
the present day number (local space density) of PCEBs is 
6.1 X 10^ (1.2 X IQ-^ PC~f ); 

ISchreiber fc Gansickd (|2003l ) estimated from the ob- 
served sample of PCEBs that 6 x 10"® < £i/pc"^ < 3 x lO"''. 
An IMRD of n(qi) oc gives the best agreement 

with this observed estimate, where we obtain 6.0 x 10"'' < 
qIvcT'^ < 5.6 X 10"^ 

Our models with qce = 0.1 also provide small space 
densities (6.0 x 10"'^ < £i/pc"'' < 1.7 x 10"^, depending on 
the IMRD) , which are similar to the observed ones. However, 
we find that in these models no PCEBs with Mwd/Mq < 
0.4 form. Low mass helium white dwarfs form via a case 
B CE phase, when the primary fills its Roche lobe on the 
first giant branch, typically when Ri ~ 15 to 100 Rq. This 
means that the orbital period at the onset of the CE phase 
is between approximately 40 and 250 d. As a consequence 
of such a low ejection efficiency, the binary cannot spiral 
in sufficiently to eject the envelope before a merger occurs. 
The lack of PCEBs with low mass white dwarfs in models 
with Qce = 0.1 is in conflict with observations; for example 
LM Com and 0137-3457 have 95 per cent and 61 per cent 
probabilities respectively of having A/wd/Mq < 0.4. 



3.5 IK Peg: A Clue to the CE Mechanism? 

As discussed above, the bulk of observed systems which 
have Af2/M0<l.O and Porb/d<l are consistent with our 
reference model A. IK Peg, on the other hand, cannot be 
explained by models with qce < 1.0, for A = 0.5. Fig- 
ure [7] shows the theoretical PCEB population with 1.1 < 
A/wd/Mq < 1.44 for a range of models, indicated in the 
bottom right hand corner of each panel. The location of IK 
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Figure 7. Calculated PCEB populations with 1.1 < J\/wd/Mq < 1.44 for a range of models indicated in the lower right hand corner 
of each plot. Here n{q{) = 1. The observed location of IK Peg is indicated by a diamond. The grey-scale bars above each plot indicates 
the number of systems per bin area. 
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Table 4. The formation rates, present day numbers and local 
space densities of PCEBs for different treatments of the CE phase, 
and for different initial mass distributions of the secondary star. 
The space densities were calculated by dividing the present day 
numbers by the Galactic volume of 5 X 10^^ pc^. 



Model 


Formation rate 


Number of systems 


Local space 




/yr-i 




density /pc— 3 


n((}i) oc gr"'™, < (?i < 1 


CEOl 


2.1 X lO""* 


3.0 X 10'' 


6.0 X 10- ' 


CE06 


1.5 X 10-3 


4.4 X IQS 


8.8 X 10-" 


A 


2.0 X 10-3 


7.3 X lOf* 


1.5 X 10-^5 


PL05 


1.6 X 10-3 


4.8 X 10'' 


9.6 X 10-" 


PLl 


1.3 X 10-3 


3.6 X lO'' 


7.2 X 10-" 


PL2 


1.1 X 10-3 


2.7 X 10'' 


5.4 X 10-" 


CT0375 


1.8 X 10-3 


6.6 X 10" 


1.3 X 10-" 


CT075 


1.7 X 10-3 


5.7 X lO'' 


1.1 X 10-5 


CT15 


1.3 X 10-3 


4.2 X 10" 


8.5 X 10-" 


DTg 


1.9 X 10-3 


6.8 X 10" 


1.3 X 10-" 


DTb 


2.5 X 10-3 


1.2 X lO'^ 


2.4 X 10-5 


nl5 


5.4 X 10-3 


2.8 X 10'' 


5.6 X 10-" 


"(Si) = 1, < gi < 1 (our reference IMRD) 


CEOl 


9.8 X 10-3 


1.5 X 10^ 


3.0 X 10-5 


CE06 


5.5 X 10-2 


1.5 X 10** 


3.0 X 10-'' 


A 


7.2 X 10-2 


2.2 X 10** 


4.4 X 10-" 


PL05 


6.8 X 10-2 


1.9 X 10** 


3.8 X 10-* 


PLl 


6.5 X 10-2 


1.6 X 10** 


3.2 X 10-* 


PL2 


6.1 X 10-2 


1.4 X 10** 


2.8 X 10-* 


CT0375 


7.0 X 10-2 


2.1 X 10** 


4.3 X 10-* 


CT075 


6.7 X 10-2 


2.0 X 10** 


4.0 X 10-* 


CT15 


6.0 X 10-2 


1.7 X 10** 


3.5 X 10-* 


DTg 


7.0 X 10-2 


2.0 X 10* 


4.0 X 10-* 


DTb 


8.7 X 10-2 


3.1 X 10** 


6.2 X 10-* 


nl5 


3.6 X 10-1 


6.1 X 10* 


1.2 X 10-3 




"(<?i) 


oc iji, < (ji < 1 




CEOl 


1.1 X 10-2 


1.9 X 10 '' 


3.8 X 10-" 


CE06 


5.8 X 10-2 


1.4 X 10* 


2.8 X 10-* 


A 


7.5 X 10-2 


2.0 X 10* 


4.0 X 10-* 


PL05 


7.7 X 10-2 


1.8 X 10* 


3.6 X 10-* 


PLl 


7.8 X 10-2 


1.8 X 10* 


3.6 X 10-* 


PL2 
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3.2 X 10-* 
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7.3 X 10-2 
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1.1 X 10" 
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Figure 8. Same as for Fig.[2]but now showing the theoretical and 
observed PCEB population with AfwD /M© < 0.4 for model PL2 
(top panel) and for model CT15 (bottom panel); in both cases 
n(gi) = 1. The grey-scale bar above the plots indicates the number 
of systems per bin area, while the bar on the right indicates the 
weightings of the observed systems. 



Peg on the M2 — log Porb plane is indicated as the diamond 
in each case. 

In order to explain the location of IK Peg, we require 
acE > 3, if A = 0.5. Model PL2 achieves acE ~ 2.9 for 
this system and, as shown in the panel labelled 'PL2' in 
Figure [7] can account for the location of IK Peg. However, 
this model also generates a low-mass cut-oflt at A/2 ~ 0.4 
M© in the PCEB population. This is a feature which is not 
consistent with observed systems, as highlighted in the top 
panel of Figure [S] LM Com and 0137-3457 have a 95 per 
cent and a 61 per cent probability of sitting in this panel. 
Thus it appears this model cannot adequately describe the 
CE phase in its present form. 

Note that models CT0375 to CT15 cannot account for 
the location of IK Peg either. From equation (O acE 1 
as M2 — > oo- As with the model PL2, the cut-off (by de- 
sign of the model) in the PCEB population at M2 = 0.15 
M0 doesn't appear to be supported by the observed sam- 
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pie. The bottom panel of Figure [8] compares the theoretical 
PCEB population with 0.4 < Mwd/M© < 0.5, with the 
corresponding observed WD+MS systems. Indeed, model 
CT15 cannot account for the observed location of, for ex- 
ample, HR Cam or 0137-3457, which have a 84 and 61 per 
cent probability of occupying this panel. 

Model DTb, where the thermal energy of the giant's 
envelope contributes to the ejection of the CE, can also ac- 
count for the location of IK Peg, while DTg cannot (see 
panels labelled 'DTb' and DTg' in Figure [Tjl. Model DTg 
shows a slight increase in the orbital period of the upper 
boundary in the theoretical PCEB population, compared to 
that shown in the panel labelled 'A', but this is not enough 
to account for IK Peg. A possible primary progenitor of IK 
Peg would fill its Roche lobe with a mass of 6 M© and a ra- 
dius of 725 R©. From Fig. [1] (dashed line), we find Ag f» 0.6. 
However, to account for the location of IK Peg, we require 
Ag « 2.2 (with acE = 1). 

Finally, the location of IK Peg can be accounted for if 
we consider an angular momentum, rather than an energy, 
budget, as shown by the panel 'nl5'. Of further interest is 
that, in contrast to the other models considered in Figure [O] 
model nl5 predicts no PCEBs with M2 < 1.0 Mq at Porb < 1 
d. Furthermore, while the PCEB population in models A to 
DTb peaks at Porb ~ 1 to 10 d, the model nl5 popula- 
tion peaks at Porb ~ 1000 d. This confirms the result by 
iMaxted et all l|2007l ). 



3.6 The Initial Secondary mass Distribution 

We now consider the impact of the initial mass ratio distri- 
bution on the theoretical PCEB populations. For our stan- 
dard model A we calculated the PCEB population for each 
initial distribution of the secondary mass, which we com- 
pared to the observed WD-I-MS systems. Figure [9] illustrates 
the differences between the resulting PCEB populations in 
the range 0.4 < Mwd/Mq < 0.5 as an example. 

For the cases n{q{) oc qi~°'^^ and IMFM2 the bulk of the 
PCEB population Ues at A/2 < 0.4 M© and Porb < 1 d. This 
is consistent with the location of the observed WD-(-MS sys- 
tems in Figure[9] The theoretical distribution with n{qi) = 1 
suggests a peak at M2 « 1 M0 and Porb 1 d, yet no sys- 
tems are observed in this region of Af2 — log Porb space. How- 
ever, this is likely to be a selection effect; PCEB candidates 
are selected according to their blue colour due to the optical 
emission from the white dwarf and/or high radial velocity 
variations. The flux from early-type secondaries will domi- 
nate over that of the white dwarf, an d henc e these systems 
go undetected. ISchreiber fc Gansickd (|2003l ) predicted that 
there is a large as yet undetected population of old PCEBs, 
wit h cool white dwarfs and lon g orbit al period. 

iRebassa-Mansergas et ahl (l2008t ) performed Monte- 
Carlo simulations to calculate the detection probability of 
PCEBs with Porb /d < 10, based on the measuremen t accu - 
racies of the Very Large Telescope (jSchreiber et al.ll2008l ). 
the SDSS, and for 1, 2 and 3a significance of the radial ve- 
locity variations. They found that approximately 6 out of 
their sample of 9 PCEBs should have Porb > 1 d, yet this 
is not the case; all of their P CEBs have Porb < 1 d, in con - 
trast with the predictions of ISchreiber fc Gansickd (|2003l ). 
Thus it is possible that the sharp decline in the population 
of PCEBs with Porb > 1 d is a characteristic of the intrinsic 



PCEB population, making the models with n{qi) oc qi~'^'^^ 
or IMFM2 more attractive. Note, however, that the local 
space density calculated for an IMRD of n{qi) (X is 
in good agreement with the observationally determined one, 
while IMFM2 is not. 

To determine if the intrinsic PCEB population does 
sharply decline for Porb > 1 d, we compared the orbital 
period distribution of our observed sample of PCEBs with 
our calculated distribution, with n{qi) = qf'^'^^, model A. 
These distributions are shown as the hashed histogram and 
the red line in the top panel, left column of Figure 1101 The 
corresponding normalised cumulative distribution functions 
(CDFs) are shown in the top panel, right column of Figure 
IIOI with the scale indicated on the right axis. Note that, in 
contrast for the observed distribution, the number of PCEBs 
in the intrinsic population gradually declines for Porb ^ 1 d, 
as opposed to a sharp decline. 

We supplemented our calculated PCEB orbital period 
distribution with the det e ction probabilities calculated by 
iRebassa-Mansergas et al.l l|2008l ') (see their Figure 7). We 
considered the detection probabilities of WD-I-MS systems 
showing radial velocity variations in their spectra with a 3a 
signifi cance ( the criterion used by Rebassa-Mansergas et al.l 
(j2003) and IRebassa-Mansergas et all (120081 ) to identify 
PCEB candidates), as detected by the SDSS (b ottom 
curve of Fig. 7 in Reb assa-Manserg as et al.l (|2008l )). As 
IRebassa-Mansergas et al.i ((2008) only calculate the detec- 
tion probabilities, V, for Porb ^ 10 d, we extrapolated 
the curve up to Porb = 100 d using the curve V = 
0.43(Porb/d)-°-3^ 

The corresponding orbital period distribution is shown 
as the green histogram in the top panels in the left and 
right columns of Figure 1101 Note that the inclusion of the 
PCEB detection probabilities has a marginal effect. Indeed, 
we still predict a gradual decline in the number of PCEBs 
with Porb > 1 d. 

To determine the likelihood that the observed and cal- 
culated PCEB orbital period distributions are drawn from 
the same parent distribution, we calculate the Kolmogorov- 
Smirnov statistic from the normalised CDF distributions, 
and therefore the corresponding significance level, ctks. A 
very small value of ctks shows that the two distributions are 
significantly different, while ctks ~ 1 shows that the two 
distributions are in good agreement. We find ctks ~ 0.11 
when comparing the observed and calculated intrinsic (red) 
orbital period distributions in the top panel, left column of 
Figure [TDj On the other hand, we find ctks = 0.35 between 
the observed and the calculated (green) orbital period dis- 
tribution with the detection probabilities included. 

We also consider the selection bias towards late-type 
secondaries by only considering the observed and theoretical 
orbital period distribution of PCEBs which have M2 /M0 < 
0.5 (middle panels of Fig. [TO]) and Af2/MQ < 0.35 (bottom 
panels). There is a better agreement between the location 
of the peaks in the observed and theoretical PCEB orbital 
period distributions. However, the theoretical distributions 
(with and without the inclusion of PCEB detection prob- 
abilities), still predict a gradual decline in the number of 
PCEBs with Porb ^ 1 d, while there is a sharp decline in the 
observed distribution. For the population of PCEBs with 
M2 < 0.5 M0 and Mq < 0.35, we find ctks = 5.6 x 10"^ 
between the observed and calculated orbital period distri- 
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Figure 9. Same as Fig. [2] but showing the population of PCEBs with 0.4 < A/wd/M© < 0.5 for different initial 
secondary mass distributions, as indicated in each panel. The grey-scale bar above each panel gives the number of 
systems per bin area, while the colour bar at the top of the page gives the weighting of each system. 



butions, with and without the inclusion of PCEB detection 
probabilities. 

Thus, we cannot reproduce the observed sharp decline 
in the number of PCEBs with Pa^h ^ 1 d, even if we take 
into account in our calculations the selection biases towards 
PCEBs with late-type secondaries, and the biases against 
the detection of PCEBs with long orbital periods. It is still 
unclear whether this sharp decline is indeed a characteristic 
of the intrinsic PCEB population, or if it is a result of further 
selection effects which have yet to be considered. 



4 DISCUSSION 

4.1 Constraining the CE phase 

We have shown that the majority of observed PCEBs (con- 
taining either a sub-dwarf or white dwarf primary) can be re- 
produced by canonical models with a constant, global value 
of acE > 0.1 for the CE ejection efflciency. The systems 
V651 Mon, FF Aqr and V1379 Aql are likely to have formed 
from a thermally unstable RLOF phase. This is contrary to 



iNelemans fc ToutI (|2005l ) who assumed them to be PCEBs, 
and attempte d reconstruct their values of qce- For the case 
of V1379 Aal lNelemans fc ToutI ([2005) could not find a so- 
lution for acE, and hence took this system as evidence for 
their '7- algorithm'. We have shown that this system could 
have formed from a thermally unstable RLOF phase. 

There is only one system, IK Peg, that is both likely 
to be a PCEB and at the same time inconsistent with the 
standard energy budget CE model. Unlike the vast majority 
of the observed sample of PCEBs, this system contains an 
early-type secondary star, and this may provide a clue to 
the ejection mechanism during the CE phase. 

Formally, the observed configuration of IK Peg requires 
Qce ^ 3. This means that a source of energy other than grav- 
itational potential energy is exploited for the ejection of the 
CE. It has been suggested that this is the thermal and ioniza- 
tion energy of the giant's envelope (Han et al. 1995, Han et 
al. 1995, Dewi & Tauris 2000, Webbink 2007). We find that 
by considering this extra energy source in our models we can 
indeed account for the location of IK P eg. However, t his is 
a concept which has be e n challenged bylHarpaj l)l998h and 
ISoker fc Harpaj l|2003h . iHarpaj (|l998l ) argue that, during 
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Figure 10. Left column: the normalised PDF of the orbital 
period distribution of the observed sample of PCEBs (hashed 
histogram) compared with the calculated distribution of the in- 
trinsic population for n{q^) = q^^'^'^ (red line), model A. The 
green line is similar to the red line, except here we also take 
into account the PCEB detect ion probability, as calculated by 
iRebassa-Mansereas et al. I ||2008|) . assuming a measurement accu- 
racy of 15 km s~^ (appropriate for SDSS spectra), to detect 3o" 
radial velocity variations. Top: all PCEBs; middle: PCEBs with 
M2 < 0.5 M0; bottom: PCEBs with M2 < 0.35 M©. Right col- 
umn: normalised CDFs of the corresponding orbital period dis- 
tributions of the observed (solid black line), with the calculated 
distributions (red and green lines), with the scale indicated on 
the right axes. Also shown are the Kolmogorov-Smirnov signifi- 
cance levels, (TKSi that the observed and calculated distributions 
are drawn from the same parent population. 



the planetary nebula phase, the opacity of the giant's en- 
velope decreases during recombination. Hence the envelope 
becomes transparent to its own radiation. The radiation will 
therefore freely escape rather than push against the material 
to eject it. 

Previous population synthesis studies have considered 
constant, gl obal values of A = 0.5 ( e.g. deKool 92, Willems & 
Kolb 2004). lDewi fc TaurisI (j2000r ) suggested that this may 
lead to an overestimation of the binding energy of the gi- 
ant's envelop e, and hence underesti mate the final PCEB or- 
bital period. iDewi fc Tauri j (|2000l ) calculated values of Ag 
for range of masses and radii, which we incorporated into 
our population synthesis code. This , however, cannot ac - 
count for IK Peg either. We note that lDewi fc TaurisI (|2000l ') 
do not calculate values of Ag for R > 600 R©. We concede 
that we may still underestimate the value of Ag for the case 
of IK Peg due to our adopted extrapolation to larger radii. 
If we do linearly extrapolate Ag for a progenitor primary of 
IK Peg with Ml = 6 Mq and Ri = 725 R© we find that 
Ag « 1.5, which is still not large enough to account for IK 
Peg. It would therefore be beneficial to calculate A in the 
mass and radius regime for PCEB progenitors. 



Rather than con s ider the C E phase in terms of en - 
ergy, iNelemans e\~sl\ (|2000l ) and iNelemans fc ToutI (|2005l ) 
describe the CE phase in terms of the angular momentum 
of the binary . Indeed, this is also a prescription favoured by 
ISokerl (I2OO4I ). Our model nl5 can account for the location 
of IK Peg. However. iMaxted et al] (|2007l ) found that for the 
7-prescription, the number of PCEBs increases with increas- 
ing orbital periods larger than approximately 100 d. Indeed, 
we also find increasing number of PCEBs with high mass 
white dwarfs with increasing orbital period, with the bulk 
of the population lying at approximatel y 1000 d. This is con - 
trast with observations by Rebassa-M ansergas et all ()2008l l. 
who found a sharp decline in the number of PCEBs with 
increasing orbital period beyond 1 d (however, see Section 
4.2). 



4.2 Observing PCEBs 

Even though we have critically examined a variety of treat- 
ments for the CE phase by comparing our models to the 
observed sample of PCEBs, we have been unable to sig- 
nificantly constrain the underlying physics. We believe this 
is mainly due to the selection effects still pervading the ob- 
served sample of PCEBs. The large majority of our observed 
PCEB sample contain late type secondaries, typically M3 
to M5. This is a consequence of the fact that until recently 
PCEB candidates were identified in blue colour surveys, such 
as the Palomar Green survey. As a result systems containing 
secondaries with early spectral types will be missed, as their 
optical flux will dominate over that from the white dwarf. 
A few exceptions include systems identified by their large 
proper motions (e.g. RR Cae) or spectroscopic binaries (e.g. 
V471 Tau). IK Peg was detected due to the emission of soft 
X-rays from the young white dwarf, which has an effective 
temperature of 40 000 K. 

The present sample of PCEBs is therefore covering 
an insufficient range in secondary masses. However, mat- 
ters are improving with the advent of the SDSS, which 
probes a large ugriz colour space. This will allow an ex- 
tra 30 PCEBs to be supplemented to the currently known 
46 systems in the foreseeable future (Gansicke 2008, private 
communication) . A complementary program is currently un- 
derway to target those WD+MS systems with cool white 
dwarfs and/or early- type secondaries in order to compen- 
sate f or the bias against such systems in the previous sur- 
veys (jSchreiber. Nebot Gomez- Moran fc Schwopell2007l ). It 
is therefore feasible that we will be able to further constrain 
our models in the near future. 

The observed sample of PCEBs also have Porb ^ 1 d, 
which can be argued to be a further selection effect; PCEBs 
are also detected due to radial velocity vari a tions of their 
spectra. However, iRebassa-Mansergas et all (j2008l ) found 
that this long-period cut-off may be a characteristic of the 
intrinsic PCEB population, rather than a selection effect. 
More precisely, we have shown in Section 3.7 that this may 
be a feature intrinsic to the population of PCEBs with late- 
type secondaries. 

We also find that an IMRD distribution of n(gi) oc 
^.-0.99 reprodu ce the local space density in ferred from 
observed PCEBs bv lSchreiber fc Gansick3 (j2003l ). as well as 
accounting for the location of the currently known sample of 
PCEBs in Ad2 — log Porb space. However, the more generally 
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preferred binary star IMRD is n(qi) = 1 (e.g. Duquennoy & 
Mayor 1991, Mazeh et al. 1992, Goldberg, Mazeh & Latham 
2003). 



5 CONCLUSIONS 

By applying population synthesis techniques we have cal- 
culated the present day population of post-common enve- 
lope binaries (PCEBs) for a range of models describing the 
common envelope (CE) phase and for different assumptions 
about the initial mass ratio distribution. We have then com- 
pared these models to the currently known sample of PCEBs 
in the three-dimensional configuration space made up of the 
two component masses and the orbital period. 

We find that the canonical model of a constant, global 
value of acE > 0.1 can account for the observed PCEB sys- 
tems with late-type secondaries. However, this cannot ex- 
plain IK Peg, which has an early-type secondary star. IK 
Peg can be accounted for if we assume that the thermal and 
ionization energy of the giant primary's envelope, as well 
as the binary's orbital energy, can unbind the CE from the 
system. IK Peg can also be explained by describing the CE 
phase in terms of the binary's angul ar momentum, accord - 
ing to the 7-prescrip tion proposed bv lNelemans et al] (|2000l ) 
and iNelemans fcTout ( 2005 ) . 

We find that the present day number (local space 
density) of PCEBs in the Galaxy ranges from 3.0 x 10^ 
(6.0 X 10"^ pc"^) for model CEOl, n(gi) oc to 
1.1 X 10^ (2.2 X 10"^ pc"^) for model nl5 and for IMFM2. 

We also find that an initial mass ratio distribution 
(IMRD) of n(gi) oc gi""'^^ gives local space densities in the 
range 6.0 x lO"'^ < fj/pc~"^ < 5.6 x 10~^, in good agreement 
with the observationally determined local space density of 
6.0 x lO"*^ < ^./pc"^ < 3.0 X 10"^ This form of the IMRD 
also predicts a decline in the population of PCEBs with 
late-type (M2 < 0.35 Mp;) secondaries, wh ich is what is ob- 
served by iR ebassa-Mansergas et al.l (|2008l ). However, while 
observations show a sharp decline in the number of PCEBs 
with orbital periods larger than 1 d, our theoretical calcula- 
tions instead predict a gradual decline. We cannot reproduce 
this sharp decline even if we take into account observational 
biases towards PCEBs with late spectral-type secondaries, 
and the selection biases against PCEBs with orbital periods 
greater than 1 d. 

Our work highlights that selection biases need to be 
overcome, especially for detecting PCEBs with early-type 
secondaries and/or cool white dwarfs. This would greatly 
advance our understanding of the CE phase. 
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